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The Constancy of the Constants of Nature: Updates 

Takeshi Chiba 

Department of Physics, College of Humanities and Sciences, Nihon University, 

Tokyo 156-8550, Japan 

The current observational and experimental constraints on the time variation of the 
constants of nature (the fine structure constant a, the gravitational constant G and the 
proton-electron mass ratio /i = mp/m^) are reviewed. 

§1. Introduction 

o \ 

CN , The title of this paper is oxymoron. We mean the title as the possibility of 

O ' variation of the quantities which are supposed to be constants. 

2^ ■ There are two aspects of the constants of nature: the system of units and the laws 

of nature. On the one hand, when we attempt to describe natural phenomena using 
I physics laws, physical quantities are expressed using physical constants: the speed 

■ of light c, Planck (Dirac) constant h{h), the gravitational constant G, Boltzmann 

constant k, electron (proton) mass me{mp), for example. When combined these 
qh' constants with the electric constant eo, we can determine all the units in the SI 

units. Thus the physical constants are closely related with the system of units. 
bJO! On the other hand, there are four fundamental forces in nature: electromagnetic, 

weak, strong, gravitational. All the phenomena in nature are described by these 
Tlj- ■ four forces. The coupling constants which describe the strength of these forces 

^ I are physical constants of fundamental importance: the fine structure constant a = 

/A-Keohc ~ 7.30 x 10^^ ~ 1/137, the Fermi coupling constant Gp/ihcf' ~ 1.17 x 
O ' 10~^GeV^^, the strong coupling constant as ~ 0.119, the gravitational constant 

\ G ~ 6.67 X 10~^^m^kg~^s~^ (a dimensionless gravitational fine structure constant 

ac = Gmp/hc ~ 5.10 x 10~^^ may also be used). Therefore, if these constants are 
not constant, then the correspondence between the experimental results and theories 
would depend on when and where the measurements are performed, which would 
result in the violation of the universality of the laws of nature. Testing the constancy 
of the physical constants thus is of fundamental importance. 
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H ' l-l- Large Number Hypothesis 

Dirac appears to have been the first who argued for the possibility of time 
variation of the constants of nature.^ As is well-known, dimensionless numbers 
involving G are huge (or minuscule). For example, the ratio of the electrostatic force 
to the gravitational force between an electron and a proton is 

e2 

Ni = - ~2xl03^ (1-1) 

GrUpnie 

where e is the electric charge, nip is the proton mass and me is the electron mass. 
Similarly, the ratio of the Hubble horizon radius of the Universe, H^^ to the classical 
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radius of an electron is 

N2= ^ c.3xW'^h-\ (1-2) 

where h is the Hubble parameter in units of lOOkms^^Mpc^^. Curiously, the two 
nearly coincide, which motivated Dirac to postulate the so-called the large number 
hypothesis.^' In his article entitled "A new basis for cosmology", he described 
Any two of the very large dimensionless numbers occurring in Nature are 
connected by a simple mathematical relation, in which the coefficients are 
of the order of magnitude unity. 
Thus if the (almost) equality A'^i = 0(1) x N2 holds always, then G must decrease 
with time G oc or the fine structure constant, a, must increase with time 

a oc t^/^Ej) since H cxt'^. 

Nowadays we know that such a huge dimensionless number like A'^i is related to 
the gauge hierarchy problem. In fact, the gauge couplings are running (however, only 
logarithmically) as the energy grows, and all the gauge couplings are believed to unify 
at the fundamental energy scale (probably string scale). The fact that A^i nearly 
coincides with N2 may be just accidental, and pursuing the relation between them 
is numerological speculation (or requires anthropic arguments). However, Pandora's 
box was opened. In the following, we mention several motivations for considering 
the variation of the constants of nature. 

1.2. Newton, Einstein, String 

Space and time in Newtonian mechanics are rigid and immutable: the absolute 
space and time which define the absolute inertial frame and exist forever even without 
matter. 

The concept of space and time in general relativity is different. The structure 
of space and time is affected by the presence of matter and thus becomes soft and 
malleable. However, the laws of physics are kept rigid: the equivalence principle 
fixes locally the laws of physics. 

On the other hand, string theory can be viewed as a framework for softening 
the laws of physics.® In string theory, the coupling constants are determined by the 
vacuum expectation values of some scalar fields and thus they are no longer constant 
at all. The situation is summarized in Table HI 





spacetime 


laws of physics 


Newton 


rigid 


rigid 


Einstein 


soft 


rigid 


String Theory 


soft 


soft 



Table I. Softening of Spacetime and Laws of Physics 

String theory is the most promising approach to unify all fundamental forces in 
nature. It is believed that in string theory all the coupling constants and parameters 
(except the string tension) in nature are derived quantities and are determined by 
the vacuum expectation values of the dilaton and moduli. However, no compelling 
mechanism how and when to fix the dilaton/moduli is known. 
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On the other hand, we know that the Universe is expanding. Then it is no 
wonder to imagine the possibihty of the time variation of the constants of nature 
during the evolution of the Universe. 

In fact, it is argued that the effective potentials of dilaton or moduli induced by 
nonperturbative effects may exhibit runaway structure; they asymptote zero for the 
weak coupling limit where dilaton becomes minus infinity or internal radius becomes 
infinity and symmetries are restored in the limit .l^'l^ Thus it is expected that as 
these fields vary, the natural "constants" may change in time and moreover the 
violation of the weak equivalence principle may be induced Ej) ,[7| (see also [8]) ,[9l) for 
earlier discussion). Moreover, the present cosmic acceleration may be induced by 
a slowly rolling light scalar field (called quintessence). Quintessence can couple to 
electromagnetic field^''' and/or gravitational field^ directly unless such couplings 
are forbidden by some symmetries. The couplings could induce the time variation 
of a and/or G. 

Hence, any detection or nondetection of such variations at various cosmological 
epochs could provide useful information about the nature of dilaton/moduli fixing 
and the coupling of the quintessence field. 

1.3. Importance as Null Tests 

We should emphasize another important aspect of checking the constancy of 
the fundamental constants: a null test. It is of fundamental importance to check 
to what extent the gravitational force obeys the inverse square law and to what 
extent the equivalence principle (the universality of free-fall) holds. Likewise, it is 
of fundamental importance to check the constancy of the fundamental constants to 
the ultimate precision. By comparing the experimental values at various epochs and 
positions, we could confirm the internal consistency of the foundation of the laws of 
physics. 

1.4. Use of Cosmology 

Cosmological observations have played important roles in testing the constancy 
of the fundamental constants, which may be evident by writing the time derivative 
in terms of a difference: 

Aa 



aAt 



(1-3) 



Therefore, in order to place a strong constraint on the time variability, one needs to 
(1) measure the constant accurately (thereby minimizing Aa/a) or to (2) measure 
for a long time (larger At). Laboratory precision experiments correspond to the 
former [Aa/a <^ 1 but At ~ 0{1) yr), while cosmological observations the latter ( 
At as much as 137 Gyr but Aa/a ~ 0{1)). 

1.5. Plan of the Paper 

In this article, we review the current experimental (laboratory, astrophysical and 
geophysical) constraints on the time variation of the constants of nature. In partic- 
ular, we consider a (sec. 2), G (sec. 3), mp/rue (sec. 4) and A (sec. 5), extending and 
updating our previous review .12]. More than ten years have passed since our previ- 
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ous review, and significant progress has been made in the experimental constraints 
on the variation (in particular thanks to the release of the WMAP data) , so it is very 
timely to update our review. See also 15) for recent reviews. For earlier expositions, 
see [T3ll .ll4 p for example. We sometimes use the units of h = c = 1 and assume 
Hq = 100/ikm/s/Mpc with h = 0.71 for the Hubble parameter and Qm = 0.27 and 
i7yl = 0.73 for the cosmological parameters taken from WMAP results ri6|) . 

§2. a 

In this section, we review the experimental constraints on the time variation of 
the fine structure constant. The results are summarized in Table HIl 
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Table II. Summary of the experimental bounds on the time variation of the fine structure constant. 
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2.1. Earth and a: Oklo Natural Reactor and Meteorites 
2.1.1. Oklo Natural Reactor. 

In 1972, the French CEA (Commissariat a I'Energie Atomique) discovered an- 
cient natural nuclear reactors in the ore body of the Oklo uranium mine in Gabon, 
West Africa. It is called the Oklo phenomenon. The reactor operated about 2 bil- 
lion years ago corresponding to the redshift z ~ 0.16 for the assumed cosmology 
{h = 0.71, Qm = 0.27, Qa = 0.73). 

Shlyakhter noticed the extremely low resonance energy [E^ = 97.3meV) of the 
reaction 

1495^ + ^^ 150^^ + 7, (2-1) 

and hence the abundance of ^^^Sm (one of the nuclear fission products of ^^^C/) ob- 
served at the Oklo can be a good probe of the variability of the coupling constants.!^ 
The isotope ratio of ^^^Sm/ is 0.02 rather than 0.9 as in natural samarium due 

to the neutron flux onto ^^'^Sm during the uranium fission. The neutron-absorption 
cross section (t{E) of the reaction Eq. (j2Tp is well described by the Breit-Wigner 
formula, 

"^^^ - 2mnE{E-ErY + r^/A' ^^'^^ 

where g is the statistical factor and F = Fn + \s the total width in terms of 
the neutron and the photon widths. From an analysis of nuclear and geochemical 
data, the operating conditions of the reactor was inferred and the thermally averaged 
neutron-absorption cross section could be estimated. The result was AE,. = E^^^° — 
E^ = (-120 ~ 90)meV[ISD and AEr = 4 ± 16 meV fTO]) . On the other hand, from 
the mass formula of heavy nuclei, the change in resonance energy is related to the 
change in a through the Coulomb energy contribution 

AEr = -1.1— MeV. (2-3) 

a 

By estimating the uncertainty in the resonance energy, Shlyakhter obtained the fa- 
mous bound a/a = 10~^^yr~^. Damour and Dyson reanalyzed the data by carefully 
estimating the uncertainty and obtained a/a = (—6.7 ~ 5.0) x 10~^^vr~^ llSp . Using 
new samples that were carefully collected to minimize natural contamination and 
also on a careful temperature estimate of th e re actors, Fujii et al. reached a tighter 
boundEH a/a = (0.2 ± 0.8) x lO'^V"^ 



Recently, the use of the Maxwell-Boltzmann distribution for low energy neu- 
tron spectrum was questioned and it is claimed that the analysis of Oklo data by 
employing a more realistic spectrum including the 1/E tail implies a decrease in a, 
Aa/a > 4.5 x 10~^ with the significance being 6a [20]) . Full-scale computations of 



*-* They noted that data is also consistent with a non-nuU result: (—4.9 ± 0.4) x 10~^^yr~^, 
indicating an apparent evidence for the time variability. However, from the analysis of the isotope 
compositions of Gd, the consistency of the Sm and Gd results supports the null results. 

**•* Note the plus sign in front of 0.2 unlike I19|l which should be consistent with their value of 
AEr and the relation Eg. p^ . 
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the Oklo reactor using modern method of reactor physics, however, show no evidence 
for such a change: Aa/a = (-5.6 ~ 6.6) x 10"^ [H]); Aa/a = (-2.4 ~ 1.1) x 10"^ 
l22p F^ The discrepancy arises because the reactor model used in '^O) is an infinite 
medium reactor model and is found to be undercritical if the reactor is made finite. 

2.1.2. Meteorites. 

Another geophysical bound on the variation of a can be obtained from the 
determination of nuclear decay rates using meteoritic data.'^''^ The isotopes which 
are most sensitive to changes in a are typically those with lowest beta-decay Q-value, 
Ql3. The isotope with the smallest Q^{= 2.66 ± 0.02keV) value is ^^^ReES]). 

The present abundances of ^^''Re and ^^''Os are given by 

( is^Re)o = ( '''Re)i exp(-A(to - ^^)), (2-4) 
( i870s)o = ( '''Os)i + ( ^s^Re), (1 - exp(-A(to - U))) , (2-5) 

where the subscripts and i refer to the present and the initial quantities and A is the 
time averag ed decay constant: A = X{t)dt/{to - U). ( ^^'^Re)i can be eliminated 
to give 

( i«^Os)o = ( ^^'Os), + ( i87Re)o (exp(A(to - U)) - l) , (2-6) 

which provides a linear relation (an isochron) between the present abundances (rela- 
tive to ^'^^Os) of ^'^^Re and ^^^Os. The slope of the linear curve determines A once 
the age to — U is independently determined. 

The ^'^^Re decay constant has been determined through the generation of high 
precision isochrons from material of known ages, particularly iron meteorites. Using 
the Re-Os ratios of IIIAB iron meteorites that are thought to have been formed 
in the early crystallization of asteroidal cores, Smoliar et al. found a ^^"^Re decay 
constant of A = (1.6666 it 0.009) x 10~^^yr~^ assuming that the age of the IIIA 
iron meteorit es i s 4.5578Gyr zt 0.4Myr which is identical to the Pb-Pb age of angrite 
meteorite^ H 

The beta-decay constant depends on as A oc [25]) . and if we assume 

that the variation of comes entirely from the Coulomb term, we find using the 
nuclear mass formula^ that AQp = —19MeV Aa/a. Hence, 

AX ,Aa , , 

— = -2.0 X 10^ . (2-7 

A a 

The bound on Aa/a over the age of the solar system ~ 4.6Gyr (z ~ 0.43) is thus 
obtained from the comparison of the ^^^Re meteoritic measurements of the time 
averaged ^^'^Re decay rate with the recent laboratory measurement^^ : Aa/a 



(2.5±16) X lO^'^ at z ~ 0.43|28Di!!!)| The difference between|28]) and|29]) comes from 



*•* Note that in the authors of \22} use the formula Eq. (|2-3|) with the opposite sign which is 
corrected here. 

Note that the authors of [H} recommend the value of A = (1.6666 ± 0.017) x 10~"yr"^ 
considering the systematic error associated with the spike calibration (see Ref.l6 in 1241) ). 

***) \Yg take the recommended value of A for the meteorites. Even if we assume 0.5% error for 
the meteoritic measurements, the bound only becomes Aa/a = (—2.5 ± 15) x 10~^. 
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the difference in the adopted laboratory measurements of the decay rate. We use a 
more recent measurement [27j) . It is to be noted, however, that the bound depends 
on the way of time dependence of the ^^^Re decay constantPSJ 

2.2. (Hyper) Fine splitting and a 

According to the Dirac equation, the energy levels of a hydrogen-like atom with 
atomic number Z are given by 

E ■ — 



2 mgC^Z^a^ nieC^Z^a^ SrUeC^Z^a'^ 
"'^'^ 2^^ n3(2j + 1) + 8^^ 



where 6j = J + ^ ~ y (j + ~ Z'^a'^, n is the principal quantum number and 
j is the quantum number associated with the total electron angular momentum. 
The fine structure, which arises due to the spin-orbit coupling, is the difference in 
energy between levels of different j for the same n. For example, for the hydrogen 
atom, £^(2P3/2) - E{2P^i2) ^ m^c^a^j'il ~ 4.53 x lO^^eV ~ /ic/(2.75cm). The 
energy levels are further split into doublets (hyperfine structure) by the coupling of 
the proton spin with the total electron angular momentum. For example, for the 
hydrogen atom, the hyperfine splitting for s states (n = l,j = 1/2) is given 
AEhf = ^rUeC^a^^gp — 5.89 x 10~^eV ~ /ic/(21.1cm), where gp is the proton 
gyromagnetic ratio. 

Since the fine structure levels depend on a, the wavelength spectra of cosmo- 
logically distant quasars provide a natural laboratory for investigating the time 
variability of a. Narrow lines in quasar spectra are produced by absorption of 
radiation in intervening clouds of gas, many of which are enriched with heavy el- 
ements. Because quasar spectra contain doublet absorption lines at a number of 
redshifts, it is possible to check for the time variation of a simply by looking for 
changes in the doublet separation of alkaline-type ions with one outer electron as 
a function of redshift.l2D''32' By looking at Si IV doublet, Cowie and Songaila ob- 
tained the constraint up to z ~ 3: |Z\a/a| < 3.5 x 10~^ l33|l . Also by compar- 
ing the hyperfine 21 cm HI transition with optical atomic transitions in the same 
cloud at z ~ 1.8, they obtained a bound on the fractional change in a up to red- 
shift z ~ 1.8: Aa/a = (a2=i.8 — ao)/" = (3-5 it 5.5) x 10~^, corresponding to 
a/a = (—3.3 it 5.2) x 10~^^yr~^ [33|) . Recently, by comparing the absorption by the 
HI 21 cm hyperfine transition (at z = 0.25, 0.68) with the absorption by molecular 
rotational transitions, Carilli et al. obtained a bound: |Z\a/a| < 1.7 x 10~^ I34p . 

Webb et al.'^Sl' introduced a new technique (called many-multiplet method) 
that compares the absorption wavelengths of magnesium and iron atoms in the same 
absorbing cloud, which is far more sensitive to a change in a than the alkaline-doublet 
method. They observed a number of intergalactic clouds at redshifts from 0.5 to 
1.6. For the entire sample (30 absorption systems of FeII,MgI and Mgll) they find 
Aa/a = (-1.09 ±0.36) x 10"^, deviating from zero at the 3 a. They noted that the 
deviation is dominated by measurements at z > 1, where Aa/a = (— 1.9ib0.5) x 10~^. 
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Moreover, Webb et al. [36|1 presented further evidence for the time variation of 
a by reanalyzing the previous data and including new sample of Keck/HIRES (High 
Resolution Echelle Spectrometer) absorption systems. The results indicate a smaller 
a in the past and the optical sample (72 systems of MgI,MgII,AlII,AlIII,SiII,CrII,FeII,NiII 
and Znll) shows a 4 cr deviation for 0.5 < z < 3.5: Aa/a = (-0.72 ± 0.18) x lO'^. 
They noted that the potentially significant systematic effects only make the deviation 
significant. 

The latest analysis of the third sample including 128 absorption systems for 
0.2< z < 3.7 gives Aa/a = (-0.543 ± 0.116) x 10"^ 38), and the sample is slightly 
updated to 143 absorption systems in 39) to yield Aa/a = (—0.573 ± 0.113) x 10~^ 
for 0.2 < z < 4.2. Again it is consistent with a smaller a in the past. The significance 
is now 4.7 a. 

More recently, Webb et al. analyzed the dataset from the ESO Very Large 
Telescope (VLT) and found the opposite trend: a was larger in the past.l^ Combined 
with the Keck samples, they claimed the spatial variation of adOj) : 

— = (1.10 ± 0.25) X 10"^(r/Glyr) cos 9, (2-9) 
a 

where r is the look-back time r = ct{z) and 6 is the angle between the direction of 
the measurement and the axis of best-fit dipole. 

If these observational results are correct, it would have profound implications 
for our understanding of fundamental physics. So the claim needs to be verified in- 
dependently by other observations. However, recent observations from VLT/UVES 
(Ultraviolet and Visual Echelle Spectrograph) using the same MM method (but 
only single species) have not been able to duplicate these results: for one group, 
Aa/a = (-0.06 ± 0.06) x 10"^ for Mgll/Fell systems at 0.4 < z < 2.3 Ml and 
Aa/a = (0.15 ± 0.43) x 10^^ for SilV systems at 1.59 < z < 2.92 "12), while for 
another group, Aa/a = (-0.04 ± 0.46) x 10"^ for Fell systems at z = 1.15 Ii3l) . 
Recently, in order to avoid the infiuence of spectral shifts due to ionization inhomo- 
geneities in the absorbers and non-zero offsets between different exposures, a new 
method of probing variability of a using pairs of FeH lines observed in individual 
exposures (called SIDAM, for single ion differential a measurement procedure) is 
proposed Using this method, a tighter bound is obtained: for Fell systems 

at z = 1.839, Aa/a = (5.66 ± 2.67) x 10"^ |MD and Aa/a = (-0.12 ± 1.79) x 10"^ 
for FeH at z = 1.15 146p . The data taken by another spectrograph HARPS mounted 
on VLT yield a similar bound: Aa/a = (0.05 ± 0.24) x 10~^ for the same Fell 
systems at z = 1.1508 H7|) . The analysis of FeH lines at z = 1.58 yields Aa/a = 
(-1.5 ±2.6) X 10-6 HE]). 

It is to be noted, however, that the analysis by Srianand et al. HT]) seems suffer 
from several fiaws^^' : only about half of the observations analyzed in Chand et 
al^^ have calibration spectra taken before and after the object exposure, although 
in their paper they mentioned that procedure has been followed for all the spectra. 
Moreover, the uncertainty in wavelength calibration in HT]) may not be consistent 
with the error in Aa/a According to the analysis of the fundamental noise 

limitation,^ the systematic errors in I4ip may be several times underestimated. 
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Fig. 1. The fine structure constant determined by quasar absorption lines at several redshifts. 
Open (blue) circles are the data from (35]) ,|36| ,|38]) (Keck) and crosses are from l41|) (VLT) filled 
circles are froml44 p J45 |l (VLT). The lower panel is the binned data of l35p J38 [) and the value 
of Aa/a in each bin is the weighted mean with its associated la error bar. The datum at 
2 ~ 0.14 is the Oklo bound [TSj. A curve is a linear (in scale factor) fit to the data: Aa/a = 
2.93 X 10"^ - 2.31 X 10"'^(1 - a) . A blue curve is a fit using I5IJ| instead of IT) : Aa/a = 
6.08(±2.22) X 10~^ - 4.85(±1.46) x 10-^(1 - a). 



Recent detailed re-analysis of Srianand et al. and Chand et al. confirms these 
concerns: flawed parameter estimation methods in a minimization analysi^^ 
(see however, [51]) ) and systematic errors in the UVES wavelength calibration^^' 0. 
A revised value correcting for respective error is Aa/a = (-0.44 ±0.16) x lO'^ {x^ 
minimizationjSSl' and Aa/a = (—0.17 it 0.06) x 10^^ (wavelength calibration).'^ 
The data taken from Keck and VLT are shown in Fig. [TJ A curve is a linear (in 

*' Recently, systematic errors in the absolute wavelength calibration of the optical spectrum 
of HIRES are identified in [53) with typical amplitudes being ±250m/s. However, their effects on 
Aa/a are found to be relatively smallPI from (-0.57 ±0.11) x 10"^ to (-0.61 ±0.11) x 10"^ 
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scale factor) fit to the data: Aa/a = 2.93(±1.80) x lO"'^ - 2.31(±1.01) x 10^^(1 -a) 
and the reduced is X^/d-o.f = 1.15 0- 

We have observed the same objects as Webb et al.'s group by Subaru telescope 
in August 2004 [55]) . The analysis of our observations would provide another inde- 
pendent useful information and help to clarify the situation. 

2.3. Laboratory Tests: Clock Comparison 

Laboratory experiments place constraints on the present day variation of a and 
are repeatable and systematic uncertainties can be studied by changing experimen- 
tal conditions, and hence such laboratory experiments are complementary to the 
geophysical or cosmological measurements. The laboratory constraints so far are 
based on comparisons of atomic clocks with ultrastable oscillators of different phys- 
ical makeup such as the superconducting cavity oscillator vs. the cesium hyperfine 
clock transitioii^ or the Mg fine structure transition vs. the cesium hyperfine clock 
transition.'^ In SI units, the second is defined as "the duration of 9192631770 peri- 
ods of the radiation corresponding to the transition between the two hyperfine levels 
of the ground state of the ^^^Cs atom".'^ A cesium atomic clock is the appa- 
ratus which tunes the microwave oscillator to the same frequency as the resonant 
absorption frequency of cesium (9192631770IIz). Such a clock comparison can be a 
probe of the time variation of a since a hyperfine splitting is a function of Za {Z is 
an atomic number) and is proportional to Za^^iij^ / iiB){fne/^p)RooFrei{Za) (where 
Frei{Za) is the relativistic correction factor, /xjv is the nuclear magnetic moment, 
fiB = eh/2mpC is the nuclear magneton, and R^o = a'^m^c/'iirh is the Rydberg 
constant). More than ten years ago, comparisons of rates between clocks based on 
hyperfine transitions in alkali atoms with different atomic number Z (H-maser and 
Hg"*" clocks) over 140 days yielded a bound on a: \a/a\ < 3.7 x 10~^^vr~^ 

Recently, by comparing a ^^^Hg+ optical clock ((^Si/2F = 0) - (^B^/2F = 2, mp = 0) 
electric-quadrupole transition at 282 nm) with a ^'^^Cs clock over 2 years, a much 
severer upper bound has been obtained: \a/a\ < 1.2 x lO^^^yr^^ 60). The electric- 
quadrupole transition of ^^^Hg"*" is expressed as -Roo-^Hg(a); where FHg(a) is a di- 
mensionless function of a. Most recent measurement over 6 years gives \a/a\ < 1.3 x 
lO^^^yr"^ 65) for ^^^Hg"^. Moreover, the comparison of the hyperfine frequencies of 
"^^■^Cs and ^''Rb atoms over nearly 5 years yields a/a = (0.045 it 1.6) x lO^^^yr^-*^ 
[6T]) . The comparison of the absolute 1S-2S transition in atomic hydrogen to the 
ground state of cesium combined with the results ofl60 p . l6T]) yields a constraint on a: 
a/a = (—0.9 lb 2.9) x 10~^^yr~^|62J). However, hyperfine frequencies are sensitive not 
only to a but also to a variation of the nuclear magnetic moment. Moreover, since 
a microwave distorts atoms, the improvements in microwave cesium clocks beyond 
10~^® are unlikely. With these motivations, recently, an optical electric quadrupole 
transition frequency at 436 nm in ^^^Yb"*" has been measured with a cesium atomic 
clock at two times separated by 2.8 years.l^ Combined with the data with those for 
optical transition frequencies in ^^^Hg"*" from [60|1 and in hydrogen from 62) gives 



*-* If we use |50| instead of l4ip . the corresponding fitting curve is Aa/a — 6.08(±2.22) x 10 ^ — 
4.85(±1.46)xI0~''(l-a) (blue curve in lower panel of Fig. 1.) and the reduced is /d.o.i = 1.53. 
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a/a = (-0.3±2.0) x IQ-^^yr-^^. Comparisons of Hg+ clockP^ with Yb+ and H 
yield a/a = (— 0.55ib0.95) x 10~^^yr~^ [65|). Also, comparisons of optical Sr clock^SS) 
with Hg+,E3 Yb+,E3J and iP* give a/a = (-3.3 ± 3.0) x lO^^'^yr-i [Mj). Recently, 
two optical clocks using A1+ ((^So) - (^Po) clock transition at 267 nmf^ and Hg+ 
are compared directly without a cesium atomic clock. The transition frequency 
depends both on the Rydberg constant and on a and can be expressed as -Roo-^Ai(«)- 
Prom the ratio of the two transition frequencies, a/a = (—1.6 it 2.3) x 10~^^yr~^ is 
obtainedjISSl' being independent of the assumptions on the constancy of other con- 
stants. The frequency uncertainties of the optical clocks are currently less than 
2.3 X 10~^^. The accuracy could soon compete with the gravitational redshifts due 
to the difference in the heights of the clocks {gAh/(? ~ 10^^^(/l/i/lcm)) so that the 
optical clocks could be used to map the gravitational potential of the earth and to 
test gravitational physics.^*^''^^' 

More recently, following the proposal of [70|1 . it is demonstrated that, instead of 
comparing atomic-clock of different atomic number, the difference of the electronic 
energies of the opposite-parity levels in two isotopes of the same atomic dysprosium 
(Dy) can be monitored directly using a radio-frequency electric-dipole transition 
between them.l^D Eight months measurements of the 3.1-MHz transition in ^^^Dy 
and the 235-MHz transition in ^^^Dy show that the frequency variation is 9.0 it 
6.7Hz/yr, — 0.6ib6.5 Hz/yr, respectively, which corresponds to a/a = (— 5.0ib3.7) x 
10-i5yr-i for the 3.1-MHz transition and a/ a = (-0.3±3.6) x lO-^^yr"! for the 235- 
MHz transition. The difference frequency gives finally a/a = (— 2.7ib2.6) x 10~^^yr~^ 
[7T]) . A unique aspect of this measurement is that the interpretation does not require 
comparison with different measurements to eliminate dependence on other constants. 
Current systematic uncertainties are at 1 Hz-level, but mHz- level sensitivity (la/al ~ 
lO^^^yr"^) may be feasible with this method. 

2.4. Cosmology and a: Big Bang Nucleosynthesis and Cosmic Microwave Back- 
ground 

2.4.1. Big Bang Nucleosynthesis. 

The process of the Big Bang nucleosynthesis proceeds as follows. When the 
temperature of the universe is greater than IMeV, protons and neutrons are inter- 
changed by the weak interaction. The neutron-to-proton number ratio {n/p) is given 
by the equilibrium condition: 

{n/p) = eM-Q/T), (2-10) 

where Q = 1.29MeV is the mass difference between neutron and proton. The equilib- 
rium is violated when the expansion rate of the universe H ~ y/GT^ becomes faster 
than the reaction rate of the weak interaction nav ~ G'j^T^, where Gp is the Permi 
constant. The balance of these two rates determines the freeze-out temperature, 

Tf ~ Gp.'^^^G^/^ ~ IMeV. (2-11) 

Por the temperature below IMeV, the number of neutrons decreases only due to the 
natural decay with the life time being 15 minutes. Deuterons, ^HeS and finally ^HgS 
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are produced by nuclear interactions of protons and neutrons at T ~ O.lMeV (the 
first three minutes). Ah the neutrons are incorporated into ^Hg and the abundance 
of ^He, Yp, is given by 1^ = 2{n/p)/[l + {n/p)]. Changes in Yp are induced by changes 
in Tf and Q. However, it is found that l^is most sensitive to changes in The 
a dependence of Q can be written st;^Siu3nSi 

Q ~ 1.29 - 0.76 X Aa/a MeV, (2-12) 

and a change in Yp is related to a change in a as 

AY AO Aa , , 

Comparing with the observed Yp{Yp = 0.249 ± 0.009)j^ a constraint on Aa/a 
is obtained: |Z\a/a| < 6 x 10~^ [TTj). A similar analysis yields a bound on AQ: 
-4 X 10-2 < AQ/Q < 2.7 x IQ-^ |78]), which can be translated into a bound on Aa 
via Eq. ([2031) as, -4.5 x 10"^ < ^q,/^, < g_7 x 10-2. 

2.4.2. Cosmic Microwave Background. 

Changing a changes the Thomson scattering cross section, ax = Sira'^ /Srn^, and 
also changes the differential optical depth f of photons due to Thomson scattering 
through f = XeUpUT, where Xe is the ionization fraction and Up is the number 
density of electrons. From the Saha equation, the equilibrium ionization fraction 
Xe^ is proportional to (me/r)^/^ gxp(— Q:2me/2r). Therefore, changing a alters the 
ionization history of the universe and hence affects the spectrum of cosmic microwave 
background fluctuations. 

The last scattering surface is defined by the peak of the visibility function, 
g{z) = e~'^^^^ dT I dz^ which measures the differential probability that a photon last 
scattered at redshift z. As explained in l79p . increasing a affects the visibility function 
g{z): it increases the redshift of the last scattering surface and decreases the thickness 
of the last scattering surface. This is because the equilibrium ionization fraction 
xf"*^, which is exponentially sensitive to q, is shifted to higher redshift (the effect of 
increase of f due to the increase of gt is minor) and because Xg more closely tracks 
x^^ for larger a. 

An increase in a changes the spectrum of CMB fluctuations: the peak positions 
in the spectrum shift to higher values of i (that is, a smaller angle) and the values of 
Ci (the angular power spectra of temperature anisotropics) increase.'^ The former 
effect is due to the increase of the redshift of the last scattering surface, while the 
latter is due to a larger early integrated Sachs- Wolfe effect because of an earlier 
recombination. Beyond the first peak, the diffusion damping of CMB fiuctuation due 
to the thickness of the last scattering surface becomes important.'^''^ The diffusion 
damping is caused by the random walk of CMB photons and hence the diffusion 
length A^) is given by Ad ~ \/\/Hf. The damping factor of CMB fluctuations is 
estimated as ~ exp(— A^^/A^) for a given wavelength A of the fluctuations. A large a 
shortens the diffusion length A/j and hence weakens the effect of diffusion damping 
and makes the values of increase. 
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The analysis of the first-year observations of CMB fluctuations from the WMAP (Wilkinson 
Microwave Anisotropy Probe) sateflite-^-'^ gives —0.06 < Aa/a < 0.01 at 2 a [8T|) . 
The five-year WMAP data combined with the CMB data sets by ACBAR, QUAD, 
BICEP, BOOMERanG and CBI and with the recent measurement of the Hubble 
constant by HST^^' improves the accuracy: —0.011 < Aa/a < 0.015 at 2 a |83|) . 
A recent analysis of the seven-year WMAP data combined with the matter power 
spectrum of Sloan Digital Sky Survey LEG yields Aa/a = -0.014 ± 0.014 at 2 u 
EH). 

2.4.3. Constraints at 30 < z < 1000 and z < 1. 

Several probes of variations in a after the epoch of last scatter have been pro- 
posed: 21 cm absorption of CMBPS and peak luminosity of type la supernovae 
(SNIa).'^ The former probes the variation for redshifts in the range 30 < z < 1000, 
while the latter at z < 1. 

After recombination (z ~ 1000) and before reionization (z ~ 30), hydrogen 
atoms are in ground state which is split into a singlet and a triplet state due to 
hyperfine splitting. The absorption of CMB at 21 cm hyperfine transition of the 
neutral atomic hydrogen is very sensitive to the variations in a. The Einstein A 
coefficient of the spontaneous emission of the 21 cm transition is proportional to a^^ 
and the brightness temperature signal of CMB at 21 cm is proportional to a^ l85|) . 
Future radio telescopes may give a constraint on Aa of 1% [85j) . 

A Type la supernova is considered to be a good standard candle, because its 
peak luminosity correlates with the rate of decline of the magnitude. Observations 
of Type la supernovae have been used to constrain cosmological parameters. The 
homogeneity of the peak luminosity is essentially due to the homogeneity of the 
progenitor mass, and this is primarily determined by the Chandrasekhar mass, which 
is proportional to G~^/'^. The peak luminosity also depends on the diffusion time 
of photons, which depends on a through the opacity. A decrease in opacity reduces 
the diffusion time, allowing trapped radiation to escape more rapidly, leading in 
turn to an increase in the luminosity. Decreasing a causes the opacity to decrease, 
which allows photons to escape more rapidly, thereby leading to an increase in the 
luminosity. Thus a smaller (larger) value of a would make supernovae brighter 
(fainter). The change in the absolute magnitude AM at the peak luminosity is 
related to the variation in a as AM — iAa/a\~^ i86|). Future experiments to observe 
distant SNIa like SNAP would reduce systematic errors to a magnitude of 0.02 mag, 
which corresponds to Aa/a < 2 x 10^^. This bound is significantly larger than the 
current bound by QSO for redshifts in the range 0.5 < z <2: Aa/a < 10~^. 

§3. G 

In this section, we review the experimental constraints on the time variation of 
the gravitational constant. For more detailed earlier review see I109p . The results 
are summarized in Table IIIII 
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3.1. Planetary motion and G 

If we write the effective gravitational constant G as G = Gq + Go{t — to), the 
effect of changing G is readily seen through the change in the equation of motion: 

(fx _ GMx _ GqMx Gq GqM x{t - tp) 

- — - -Q-^— -2 • y^- I 

Thus the time variation of G induces an acceleration term of secular type in addition 
to the usual Newtonian and relativistic ones, which would affect the motion of bodies, 
such as planets and binary pulsar. 

A relative distance between the Earth and Mars was accurately measured by 
taking thousands of range measurements between tracking stations of the Deep Space 
Network and Viking launders on Mars. From a least-squares fit of the parameters of 
the solar system model to the data taken from various range measurements including 
those by Viking landers to Mars (from July 1976 to July 1982), a bound on G is 
obtained: G/G = (2 ± 4) x IQ-^^y^-i 

Similarly, Lunar-Laser-Ranging measurements have been used to accurately 
determine parameters of the solar system, in particular the Earth-Moon separa- 
tion. From the analysis of the data from 1969 to 1990, a bound on G is ob- 
tained: G/G = (0.1 ± 10.4) X 10-i2yr"-^ ISSI); while from the data from 1970 to 
1994, G/G = (1 ± 8) X 10-i2yj,-i ^ Recent analysis using the data up to April 
2004 yields G/G = (4 ± 9) x lO-^^yp-i ^ ^he uncertainty for G/G is improving 
rapidly since the sensitivity for the observations depends on the square of the time 
span. 





redshift 


AG/G 


G/G(yr-i) 


Viking Lander Rangin^^ 







(2 ±4) X 10"^^ 


Lunar Laser Rangin^^ 







(4 ±9) X 10"^^ 


Double Neutron Star Binarjl^ 







(1.10 ± 1.07) X 10"" 


Pulsar-White Dwarf Binarji^ 







(-5± 18) X 10"^^ 


Helioseismologj-'^ 







< 1.6 X 10"^^ 


White Dwarf Luminosity FunctiorP^ 







< 1.8 X 10"^^ 


Neutron Star MasS^iii' 


0-3-4 




(-0.6 ±2.0) X 10"^^ 


Gravochemical Heatin^^^ 







< 4 X 10"^° 




10^ 


-0.3 ~ 0.4 


(-2.9 ~ 2.2) X 10"" 


BBN+CMdi^ 


10^ 


-0.15 ~ 0.21 


(-1.5 - 1.1) X 10"" 


BBN+CMtf^ 


10^ 


-0.10 ~ 0.13 


(-0.95 ~ 0.73) X 10"" 




10^ 


< 0.05 


< 3.6 X 10"^^ 



Table III. Summary of the experimental bounds on the time variation of the gravitational constant. 

AG/G = (Gthen — Gnow ) /G*now . 



3.2. Binary Pulsar and G 

The timing of the orbital dynamics of binary pulsars provides a new test of time 
variation of G. To the Newtonian order, the orbital period of a two-body system is 
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given by 

a3 \ 27rf3 



where a is the semi-major axis, i = r'^cf) is the angular momentum per unit mass, m 
is a Newtonian-order mass parameter, and e is the orbital eccentricity. This yields 
the orbital-period evolution rate 

Pb G i m , , 

Damour, Gibbons and Taylor showed that the appropriate phenomenological limit 
on G is obtained by 

where dPb represents whatever part of the observed orbital period derivative that is 
not otherwise explained!^ Prom the timing of the binary pulsar PSR 1913-1-16, a 
bound on G is obtained: G/G = (1.0 ± 2.3) x lO'^yr-^ EI]) (see alsoESD) However, 
only for the orbits of bodies which have negligible gravitational self-energies, the 
simplifications can be made that Pb/Pb is dominated by —2G/G term. When the 
effect of the variation in the gravitational binding energy induced by a change in 
G is taken into account, the above bound is somewhat weakened depending on 
the equation of state. This may not be concern to neutron star - white dwarf 
binaries such as PSR B1855-F09 {G/G = (-9± 18) x 10~i2yj,-i) [g^ ^nd J0437-4715 
{G/G = (-5 ± 18) X 10-i2yr"^) ESI)- 

3.3. Stars and G 

Since gravity plays an important role in the structure and evolution of a star, 
a star can be a good probe of the time variation of G. It can be shown that the 
luminosity of a star is roughly proportional to G^ if free-free transition dominates 
the opacity.!^ Increasing G is effectively the same, via the Poisson equation, as 
increasing the mass or average density of a star, which increases its average mean 
molecular weight and thus increases the luminosity of a star and hence decreases 
its lifetime. Since a more luminous star burns more hydrogen, the depth of convec- 
tion zone is affected which is determined directly from observations of solar p-mode 
(acoustic wave) spectra.'^ Helioseismology enables us to probe the structure of the 
solar interior. Comparing the p-mode oscillation spectra of varying-G solar mod- 
els with the solar p-mode frequency observations, a tight bound on G is obtained: 
\G/G\ < 1.6 X 10-i2yj.-i 98])^ 

The balance between the Fermi degeneracy pressure of a cold electron gas and 
the gravitational force determines the famous Chandrasekhar mass 

Mch ^ G-^/^mf, (3-5) 

where rrip is the proton mass, which is the upper bound of the masses of white 
dwarfs. White dwarfs are long-lived objects (~ 10 Gyr) and their inner cores are 
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almost degenerate, and hence even small variations of G can affect their structure 
and evolution. Moreover, white dwarfs do not have nuclear energy sources and their 
energy is of gravitational and thermal origin. The cooling process of white dwarfs 
is now well understood, including the energy release due to ^^Ne sedimentation in 
the liquid phase and due to C/0 phase separation on crystallization in the core.!® 
The decrease in G (larger G in the past) accelerates the cooling of white dwarfs. By 
comparing the white dwarf luminosity function measured in the open cluster NGC 
6791 with the simulated luminosity function, using the observed distance modulus 
to break the degeneracy between the age of the cluster and the effect of G, a tight 
bound on G is obtained; G/G > -1.8 x lO'i^ fTOO]) . 

Mch sets the mass scale for the late evolutionary stage of massive stars, including 
the formation of neutron stars in core collapse of supernovae, and it is thus expected 
that the average neutron mass is given by the Chandrasekhar mass. Measurements 
of neutron star masses and ages over < z < 3 ~ 4 yield a bound on G, G/G = 
(-0.6 ±2.0) X 10~^^vr"^ IT0T)) . 

Recently, a new method for constraining G is proposed using the surface tem- 
peratures of neutron stars (dubbed "gravo chemical heating" ) An increase (or 
decrease) in G induces the compression (or expansion) of the star. Since the chemical 
potentials depend on the density, the system interior to the star departs from the 
beta equilibrium state, which increases the chemical reaction rates so as to reach a 
new equilibrium state, dissipating energy as internal heating and neutrino emission. 
Comparing the ultrainfrared observation of the surface temperature of the millisec- 
ond pulsar (PSR J0437-4715), upper limits on G are obtained: \G/G\< 2xl0-iVr"^ 
if direct Urea reaction operating in the neutron star core is allowed, while \G/G\ < 
4 X lO^^'^yr"^ if only modified Urea reactions are considered .1^^ 

3.4. Cosmology and G: Big-Bang Nucleosynthesis and Cosmic Microwave Back- 
ground 

3.4.1. Big Bang Nucleosynthesis. 

The effect of changing G on the primordial light abundances (especially ^Hg) 
is already seen in Eq. (|2-10p and Eq. ()2Tip : an increase in G increases the expansion 
rate of the universe, which shifts the freeze-out to an earlier epoch and results in a 
higher abundance of ^Hg. In terms of the "speed-up factor" (the ratio of the Hubble 
parameter to that in the Standard Big Bang Nucleosynthesis), ^ = H/Hsbbn, Yp is 
well fitted bjfl® 

Yp ~ 0.244 + 0.074(^2 - 1). (3-6) 

If 1^ was between 0.22 and 0.25, then —0.32 < AG/G < 0.08, which corresponds 
to G/G = (—0.55 ~ 2.2) x 10~^^yr~^. A similar (more conservative) bound was 
obtained in I104p : -0.3 < AG/G < 0.4. Combining the determination by WMAP 
of l?^/!^ and recent measurements of primordial deuterium abundance (but with- 
out Helium and Lithium abundance), a slightly tighter constraint is obtainecP® 
: -0.15 <Z\G/G< 0.21. By combining WMAP value of Qsh"^ and recent re- 
sults of the reanalysis of helium abundance,!^ a similar bound has been obtained: 
-0.10 < AG/G < 0.13[77D. It should be noted however that the analysis by WMAP 
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team assumed the Einstein gravity and the effect of changing G is not included to 
determine Qsh'^- Hence, these analyses are not consistent and should be made in 
the context of scalar-tensor gravity (or its variants) consistently. 

3.4.2. Cosmic Microwave Background. 

Changing G changes the Hubble parameter and hence changes the size of horizon: 
H^^ oc G^^/^, which results in the change of both the location and the amplitude 
of acoustic peaks through the projection effects and the diffusion damping scale .1^^ 
For example, an increase in G shifts the peak positions in the spectrum toward higher 
values of L A larger G makes the diffusion length — Hf shorter and thus 
weakens the diffusion damping at the peak positions because the peak positions also 
depend on (A oc H~^) and hence the damping factor exp(— A|)/A^) decreases 
less. 

Recently, anisotropics in the cosmic microwave background have been measured 
up to £ < 800 by WMAP sateUite.tiS From the analysis of the WMAP data within 
the context varying G model, the variation of the gravitational constant at the 
recombination epoch is constrained as,l^^ AG/G < 0.05. 

3.5. Measuring Gq: Recent Developments 

In 1798, Cavendish carried out experiments to measure Gq by using a torsion 
balance apparatus (proposed and constructed by John Michell), which has become 
known as the Cavendish Experiment .1® The method is still used basically in mea- 
suring Gq. 

The torsion balance consists of a dumbbell suspended from the middle by a 
thin fiber. The dumbbell consists of two small masses (m) fastened to a thin rod of 
length 2i. When a large pair of masses (M) are brought into proximity to the smaller 
masses, the dumbbell rotates by an angle <po and comes to a halt. From the balance 
between the torsion of the fiber (the torsion coefficient D) and the torque due to 
the gravity force between m and M (the distance r), we have D^q = 2GoMm£/r'^ . 
When the large masses are removed from the set up, the dumbbell begins to oscillate 
because of the restoring force of the fiber. The period of the oscillation, T, is given 
by T = where /(= 2mi'^) is the moment of inertia of the dumbbell. 

Eliminating D, Gq is thus determined by G = Air'^r'^ipQi/MT'^ . 

No laboratory measurements of G/G has been performed recently (see I109p 
for older laboratory experiments). This is mainly because the measurements of 
the present gravitational constant Gq itself suffer from systematic uncertainties and 
have not been performed with good precision. For example, due to the uncertain- 
ties associated with the frequency dependence of the torsion coefficient of fibers, the 
recommended value of Gq by CODATA (Committee on Data for Science and Tech- 
nology) became worsened from G = (6.67259 it 8.5 x 10~^) x 10~^^m^kg^^s~^ in 
1986 to (6.673 ± 1.0 x 10"^) x 10"^^ in 1998.^^ 

Gundlach and Merkowitz measured Gq with a torsion-balance experiment in 
which string-twisting bias was carefully eliminated The result was a value of 
Go = (6.674215 ± 0.000092) x lO'^m^kg-^s-^. Recently, however, the measure- 
ment of G with a torsion-strip balance resulted in Gq = (6.67559 it 0.00027) x 



18 



Takeshi Chiba 




6.665 



C0D86 C0D98 00 01 020306 COD02 C0D06 CODIO 



Fig. 2. Experimental results on the measurements of G and the CODATA recommended value of 
G. 



10~^^m^kg~^s~^, which is 2 parts in 10^ higher than the result of Gundlach and 
Mer kowitz .1^^^ Probably the difference is still due to systematic errors hidden in one 
or both of the measurements. On the other hand, recent beam balance measurement 
(Go = (6.67407± 0.00022) x lO^^^m^kg-^s-^)!!!]) and torsion balance measurement 
(Go = (6.67387 ± 0.00027) x lO^^^m^kg-is^^) fm)) of G are consistent with the re- 
sult by Gundlach and Merkowitz. A new measurement of G was made using a beam 
balance.lHSJ The measured value, Go = (6.674252 ± 0.000109) x lO^^^m^kg^^s-^ is 
consistent with that by Gundlach and Merkowitz. The recommended value of G 



*' However, two recent determinations of G, one by comparing the time it took for a torsion 
pendulum to swing past masses placed at varying distances from ilP^ and another by using a laser 
interferometer to measure the displacement of pendulum bobs by various masses,'^^ give values 
significantly deviated from the value by Gundlach and Merkowitz: Go = (6.67349 ± 0.00018) x 
lO-^m^kg-^s-^ inni) and Go = (6.67259 ± 0.00085) x lO^^m^kg-^s-^ [TT7|. The source of the 
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by CODATA is being improved: G = (6.6742 ± 1.0 x 10"^) x lO^^^m^kg^^s"^ fo^ 
CODATA 2002,E3 while G = (6.67428±6.7x 10"^) x IQ-iim^kg-^s-^ for CODATA 
The values of G as a function of the year of measurements are shown in 

Fig. El 

Moreover, conceptually different experiments of measuring G using a gravity gra- 
diometer based on cold-atom interferometry were p er for med Freely falling 
samples of laser-cooled atoms are used in a gravity gradiometer to probe the field 
generated by nearby source masses. A measured value of G is, Go = (6.667 ± 0.01 1 ± 
0.003) X 10~^^m^kg~^s~'^ il22|) . It may be possible to push the measurement accuracy 
below 10"^. 

As the accuracy of the measurements improves, it may be possible to place a 
bound on the present-day variation of G. Although the current accuracy of G is more 
than 6 digits worse than the solar system experiments or cosmological constraints to 
place a constraint on the present-day G, the situation will be changed after a century 
since the accuracy improves by one digit during these ten years. It is important to 
pursue laboratory measurements of G/G since they are repeatable and hence are 
complementary to astrophysical and geophysical constraints. 

§4. Proton-Electron Mass Ratio 

In this section, we briefly mention the experimental constraints on the variation 
of the electron-proton mass ratio, /j, = nip /me. The results are summarized in Table 

4.1. Molecular lines and // 

Thompson noted the different dependence of the electronic, vibrational, and 
rotational energy levels on and first pointed out the possibility that the presence of 
cosmological evolution in can be tested by using observations of molecular hydrogen 
in quasar absorption systems.l^^ In the Born-Oppenheimer approximation, the 
molecular hydrogen levels depend on nip and can be written as 

E = E,i,, + ^ + ^. (4-1) 

Thus, the energy shift in any vibration-rotation transition j in the Lyman series has 
the form 

AEj = aeiec + bj/^ + Cj/fj,, (4-2) 
and the difference in energy between two transitions is 

AEi - AEj ~ bij/^ + Cj/^i. (4-3) 

inconsistency is currently unknown, and these new values may make the next recommended value 
of G by CODATA decrease and make the uncertainty larger. 

*' The 2010 CODATA recommended value of G is G = (6.67384 ± 8.0 x 10~*) x 
10-"m3kg-is-2.E23 
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Table IV. Summary of the experimental bounds on the time variation of the mass ratio of electron 
and proton, fi ^ rrip/me. Afj./fi = (/ithcn - Mnow)/Mnow. 



Hence, to lowest order, a change in fj, induces a change in AEi — AEj-. 

6j, ^ J{AE,-AE,) ^ 5v 2AE, 

/i AEi - ^Ej c AEi - AEj ' ^ ^ 

where 6v is the mean offset, compared to the laboratory value, of the energy difference 
between the two sets of lines, when that offset is represented as a velocity difference. 
Based on this method, Pagel first obtained a constraint on /i, < 5 x 10~^^yi~^ 

from the comparison of different redshifts determined by neutral hydrogen molecule 
and by heavy ion absorption lines. The ratio of the hyperfine 21 cm absorption 
transition of neutral hydrogen to an optical resonance transition is dependent on 
Qplia^^ where Qp is the proton gyromagnetic ratio. Comparing the measured redshifts 
of 21 cm and optical absorption a constraint on the change in /i is obtained at 
z = 1.7764 [33|) : A/j,/^ = (—0.7 ± 0.6) x 10~^ assuming that a and gp are constants 
(note the different definition of /i used there and we changed the 95% confidence 
limits given in I33p into Icr ones). Various observational constraints obtained so far 
are summarized in Table IIVI 

Recent measurements of H2 lines of Lyman and Werner bands at z = 2.5947 
and 3.0249 toward the quasars Q0405-443 and Q0347-383 with VLT/UVES indicate 
a systematic shift of ^ in the past, Afi/fj, = (1.65 ± 0.74) x 10~^ 1129]) . but the results 
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depend on the laboratory wavelengths: the above value is for wavelengths derived 
from a direct determination using laser techniques, while = (3.05ib0.75) x 10~^ 

for those derived from energy level determinationP29 In general, a measured i-th 
molecular line wavelength Aj in an absorption system at redshift Zabs is given bjff^S 

Xi = XUl + Zabs)il + K,Afi/fi), (4-5) 

where A is the laboratory transition wavelength and Ki = din Aj /d In fi is the sen- 
sitivity coefficient. With an improved calculation of sensitivity coefficients Ki and 
new accurate laboratory spectroscopic measurements, the reanalysis of the datsP^ 
strengthens the case for a larger fi in the pastPS Afi/fj, = (2.44 it 0.59) x 10~^. 
However, it is pointed out that the techniques used to calibrate the wavelength 
scale of the Ultraviolet and Visual Echelle Spectrograph (UVES) on VLT pro- 
duce calibration errors.^^ Reanalysis of the spectra by using the improved wave- 
length calibration techniques and improved fitting procedures yields a constraint 
on An/n as Aji/ji = (2.6 it 3.0) x 10~^, which is consistent with no variation 
at 1 (7 llSip . Recent observations of a number of hydrogen lines (H2 and HD) 
in the spectrum of J2123-0050 with Keck/HIRES at z = 2.059 yield a constraint 
on Afi: Afi/n = (5.6 ± 5.5(stat) ± 2.9(sys)) x 10"^ [IM])- The analysis of the 
spectrum of the same object observed with VLT/UVES gives a similar constraint: 
Afi/fi = (8.5 ± 3.6(stat) ± 2.2(sys)) x 10~^ 1133 1). The analysis of a new spectrum of 
Q0528-250 with VLT/UVES yields Afi/fi = (0.3 ± 3.2(stat) ± 1.9(sys)) x lO'^ fTMD . 

Recently, the method of using 18 cm OH lines has been proposecP^''^^''^^ to 
avoid possible systematic errors from multiple species which may have systematic 
velocity offsets. The ground "^11^/2 J = 3/2 rotation state of OH is split into two 
levels by A doubling and each of these A doubling states is further split into two 
hyperfine states. Transitions between these levels lead to four spectral lines with 
wavelength ~ 18 cm. Transitions with AF = are called the main lines, with 
frequencies of 1665.4018 and 1667.3590 MHz, while transitions with AF = 1 are 
called satellite lines, with frequencies of 1612.2310 and 1720.5299 MHz. Since the four 
OH lines arise from two very different physical processes, yl-doubling and hyperfine 
splitting, the transition frequencies have different dependences on the fundamental 
constants, a and fi and the proton gyromagnetic ratio gp. Therefore, measurements 
of these lines enable us to constrain variations in a and /i from a single species. 
The radio observations of two satellite lines at 2; = 0.247 yield AX/X = (2.2 it 
3.8) X 10~^ for X = gp^fia'^Y'^^ I138p . Assuming that the variations of a and the 
proton gyromagnetic ratio are small, a change in fi is constrained as. An/ (jl = (1.2 zt 
2.0) X 10~^ 1138]) . Deep Westerbork Synthesis Radio Telescope and Arecibo Telescope 
observations of these satellite lines yield AX/X = (—1.18 it 0.46) x 10~^, suggesting 
2.6cr evidence for a change in X I139p . The limiting cases, assuming that only a 
or 11 changes, are Aa/a = (-3.1 ± 1.2) x 10"*^ and A^/n = (-6.2 ± 2.4) x 10"^ 
I139p . All four 18 cm OH lines have recently been detected at 2; = 0.765 with low 
signal-to-noise ratio,'^^ which, when combined with HI 21 cm lines at z = 0.685, 
yields a constraint of \ Afj./fj.\ < 1.4 x 10~^ I140p . 

It is pointed out recently that the inversion transition frequencies of ammo- 
nia (NH3) are significantly sensitive to the variation of By comparing the 
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inversion spectrum of NH3 at z = 0.6847 (toward B0218+357) with rotational spec- 
tra of other molecules (CO,HCO+,HCN), a strong constraint is obtained: A^/^ = 
(—0.6 lb 1.9) X 10~^ I14ip . More detailed comparison of the NH3 inversion transi- 
tions with HCO"*" and HCN molecular rotational transitions gives a stronger con- 
straint on An: A^/n = (0.74 ± 0.47(stat) ± 0.76(sys)) x 10"'^ and a 2a constraint 
of \Afi/ fi\ < 1.8 X 10~^ I142p . The bound is recently updated further by comparing 
the NH3 inversion lines with molecular rotational lines (CS, H2CO) by using the 
Green Bank Telescope: Afi/fi = (—3.5 it 1.2) x 10~^.U^ Also, from the compari- 
son of the NH3 inversion transitions with HC3N molecular rotational transitions at 
z = 0.89 observed with the Effelsberg radio telescope, a similar bound is obtained: 
Afi/fi = (0.08±0.47) X 10-^ and a 3a constraint of \An/n\ < 1.4 x 10"^ [HI]). More- 
over, from the spectral observations of molecular cores in the disk of the Milky Way 
in molecular transitions of NH3 and HC3N at the Effelsberg radio telescope, a statis- 
tically significant velocity offset 23ib4(stat) ib3(sys)m/s between the radial velocities 
NH3 and HC3N is found.'^^SJ When interpreted in terms of the local (spatial) varia- 
tion of n, this implies a tentative signal of Afi/fj. = (—26 ± l(stat) it 3(sys)) x 10~^ 
I145p . where Afi = (^MiikyWay — IJ-iab)/ fJ-iah- However, since the number of sources is 
small and a different velocity offset is observed at a different epoch of the observa- 
tions, there may exist unaccounted-for systematic effects. 

Very recently, torsion- vibrational frequencies of methanol (CH3OH) are found to 
be far more sensitive to the variation of ' Using the published data of observing 
narrow emission lines of the methanol masers in the Milky Way, the local (spatial) 
variation of fj, is constrained as Afi/ = (—11 it 17) x 10~^ I147p . 

4.2. Laboratory Tests: Clock Comparison 

Laboratory limits on the variations of fi are also obtained by comparison of 
atomic clock^^'l^'^^ as explained in Sec. 2.3. As for the constraint using molecular 
clocks, recently, from the comparison of the frequency of a rovibrational transition in 
SFg with the hyperfine transition in Cs, combined with l65p to break the degeneracy 
with variations of a and nuclear magnetic moment, a constraint oi fi/fi = (—3.8 it 



Finally, we briefly comment on the potential variability of the cosmological con- 
stant (or dark energy) because in the runaway scenario of dilaton or moduli (p, a/a 



5.1. Evidence for 71 > 

There are two arguments for the presence of dark energy. The first indirect 
evidence comes from the sum rule in cosmology: 



5.6) X 10-^^yr-i is obtained.l^H} 



§5. A or Dark Energy 



and G/G would close to </)/</) E]). 




(5-1) 



where f2i = SirGpi/SHQ is the density parameter of the i-th energy component, pi. 
The density parameter of the curvature, fix, is defined by fix = —k/a^HQ. Since 



Constancy of the Constants of Nature 



23 



the current observational data indicate that matter density is much less than the 
critical density f2M < 1 and that the Universe is flat, we are led to conclude that 
the Universe is dominated by dark energy, f2£)E = I — f^M — > 0. 

The second evidence for dark energy is from the observational evidence for the 
accelerating universe-'^^^' PI : 



I ((2m{1 + zf + (1 + 3w)f2DE{l + zf^^+'^A > 0, (5-2 



where w is the equation of state of dark energy, w = pde/pde- Since distance 
measurements to SNIa strongly indicate the Universe is currently accelerating, the 
Universe should be dominated by dark energy with negative pressure {w < 0). We 
note that another argument for negative pressure comes from the necessity of the 
epoch of the matter domination. 

5.2. Supernova and A 

A current bound on the equation of state of dark energy from supernova data 
(580 supernovae) is Iili — 1| < 0.07 115U() . Future observations of high redshift super- 
novae/galaxies/clusters/BAO would pin down the bound on u) to |u) — 1| < 0.01. The 
extent of the time variation of dark energy density is readily seen from the equation 
of motion: 

^ = -3{l + w)H. (5-3) 

PDE 

§6. Conclusion 

A short account of the experimental constraints on the time variability of the 
constants of nature (a, G and fi) was given. Since there are some theoretical motiva- 
tions for the time variability of the constants of nature and the implications of it are 
profound, it is worth examining whether the constancy of the constants of nature is 
just a very good approximation. 

Let's keep shaking the pillars to make sure they're rigid 
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